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ABSTRACT 

b 
o 

<N ' The process of star formation in interstellar molecular clouds is believed to be 

controlled by driven supersonic magnetohydrodynamic turbulence. We suggest 
that in the inertial range such turbulence obeys the Kolmogorov law, while in the 
dissipative range it behaves as Burgers turbulence developing shock singularities. 
On the base of the She-Leveque analytical model we then predict the velocity 
power spectrum in the inertial range to be Ef. ~ fc -L74 . This result reproduces 
the observational Larson law, (uf) ~ ^.74-0.76^ [L arsori; MNRAS 194 (1981) 
^ \ 809] and agrees well with recent numerical findings by Padoan and Nordlund 

[astro-ph/0011465]. The application of the model to more general dissipative 
structures, with higher fractal dimensionality, leads to better agreement with 
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recent observational results. 

Subject headings: MHD: Turbulence — ISM: dynamics — stars: formation 

1. Introduction 



X 

It was recently argued on both observational and numerical grounds that star forming 
regions of interstellar molecular clouds are governed by super-sonic and, possibly, super- 
Alfvenic turbulence, see, e.g., (Padoan & Nordlund 1999), (Padoan & Nordlund 2000, here- 
after PN (2000)), and a review by Elmegreen (2001). The turbulence is driven on large 
scales by supernovae explosions and energy is then transfered to smaller scales via a turbu- 
lent cascade, forming a hierarchy of dense clumps. It is still unclear whether such turbulent 
fragmentation is crucial on small scales, where Jeans-unstable density cores collapse and stars 
are formed. However, it seems reasonable that at least at the initial stage of a clumpy struc- 
ture formation, turbulent fragmentation is the definitive process. This assertion, stemming 
from the work by Larson (1981), was recently confirmed in the number of high-resolution 
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numerical simulations (PN 2000; Padoan et al 2000; Klein, Fisher, & McKee 2000; Klessen 
2001a,b; Geyer & Burkert 2001; Williams 2001; Mac Low et al 2001). 

Observations suggest that the Mach number of turbulent motion, M, can be greater 
than 10, and the Alfvenic Mach number, M , can be greater than 1 [see, e.g., Klessen 
(2001a); Williams (2001)]. Until recently, supersonic turbulence (both Navier-Stokes and 
MHD) has not received proper theoretical attention. In a series of papers, Porter, Woodward, 
& Pouquet (1998) analyzed numerically decaying turbulence with initial Mach number of 
the order of 1. It has been observed in large resolution runs (up to 1024 3 ) that the spectra 
of both the compressible and incompressible parts of the velocity field approximately follow 
the Kolmogorov value, ~ k 2 \uk\ 2 ~ k~ 5 ^ 3 . However, decaying turbulence is different from 
forced turbulence in many aspects. To mention just a few, we note that a supersonic motion 
forms shocks and quickly, on a crossing time, dissipates in decaying runs, while it can be 
sustained in forced ones. Also, it has been demonstrated by Smith, Mac Low, & Zuev (2000); 
Smith, Mac Low, & Heitsch (2000) that, in a decaying case, most energy is dissipated in a 
large number of weak shocks contrary to a forced case where the largest shocks dissipate most 
of energy. In the present paper we consider supersonic, driven turbulent systems, stressing 
that they differ qualitatively from their subsonic, decaying counterparts. 

In the last two years there appeared a number of papers analyzing numerically forced 
supersonic turbulence both with and without magnetic fields. Porter et al (1999) investi- 
gated forced non- magnetized turbulence with Mach number of the order of 1, and observed 
no difference in power spectra with the unforced runs. However, when PN (2000); Padoan et 
al (2000) simulated supersonic MHD turbulence (M ~ 10, M a ~ 3), they found the velocity 
spectrum, /c _/3 , with approximate value f3 = 1.8. This spectrum is steeper than the Kol- 
mogorov one, which indicates strong intermittency effects. Correspondingly, velocity fluctu- 
ations scale with distance according to (uf) ~ I 13 ' 1 . The steeper-than-Kolmogorov spectrum 
was linked to the supersonic nature of turbulence by Larson (1979, 1981) on observational 
grounds. 

Our interest in supersonic turbulence is also motivated by the argument of PN (2000) 
that the spectral exponent, (3, may be directly related to the exponent of the mass distri- 
bution of collapsing cores, N{m) ~ m -1-5 , as 5 = 3/(4 — (3). This suggests that the initial 
mass function (IMF) could be explained from the basic properties of turbulent fragmentation, 
without tunable parameters. The fact that supersonic MHD turbulence leads to sustaining 
of shock turbulence, to shock fragmentation, and to establishing a certain universal density 
distribution has also been recently demonstrated by Boldyrev & Brandenburg (2001) in a 
one-dimensional solvable Burgers model. 

In this paper we present a theoretical model of driven supersonic turbulence, incorpo- 
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rating both Kolmogorov and Burgers pictures in different parts of the phase space. We argue 
that due to mostly solenoidal character of such turbulence, the characteristic times of energy 
cascade in the inertial interval scale as in the Kolmogorov turbulence, while the dissipative 
structures are completely different. Instead of filaments, as in an incompressible case, they 
can appear as sheets which is more consistent with Burgers turbulence. In Sec. 2 we demon- 
strate that the standard She-Leveque model (She & Leveque 1994; She & Waymire 1995), 
which links the most singular turbulent structures with turbulent spectra, has a solution 
corresponding to sheet-like dissipative structures, which reproduces the velocity power spec- 
trum with exponent (3 = 1.74, close to the observational and numerical values. In Sec. 3 we 
generalize the results to more realistic, fractal dissipative structures, which leads to better 
agreement with recent observational results. 



2. Kolmogorov-Burgers model of supersonic turbulence 

At first sight, turbulence with small pressure should behave in the same way as Burgers 
turbulence, the theory of which was substantially developed during the last few years (Polyakov 
1995; Yakhot k Chekhlov 1996; Boldyrev 1997; E et al 1997; Gotoh & Kraichnan 1998; E 
& Vanden Eijnden 1999; Verma 2000; Frisch & Beck 2000). However, this is true only in 
one- and two-dimensional cases; in a three-dimensional case, the behavior of a compress- 
ible fluid is qualitatively different from Burgers turbulence. The main difference is vorticity 
generation, an effect completely analogous to magnetic field generation existing in 3D and 
non-existing in 2D. Indeed, the vorticity equation, 

d t n + (u • v)n - (n • v)u + (v • u)n = uAn, (i) 

where the vorticity is fi = V x u, coincides with the induction equation for a magnetic 
field. Numerical experiments show that vorticity is generated quite effectively. In decaying 
turbulence with Mach numbers of the order of 1, simulated by Porter, Woodward, & Pouquet 
(1998), it was found that the turbulence was mostly solenoidal. If one decomposes the 
velocity field into the solenoidal part, V ■ u s = 0, and the compressible part, V x u c = 
0, their ratio was observed to be 7 = (u 2 c ) / {u 2 s ) ~ 0.1 in the inertial range. A pressure 
term ensuring incompressibility in subsonic turbulence, turned out to be unimportant in 
supersonic dynamics: energy transfer over scales due to the pressure term was only 3%. The 
subsequent forced runs by Porter et al (1999) revealed qualitatively the same results. In 
the case of forced turbulence with large Mach numbers (M ~ 10, M a ~ 3) and a solenoidal 
large-scale force, simulations by PN (2000) also demonstrated that in the inertial interval 
this ratio is small, 7 < 0.2, but increases towards the dissipative region. This result is 
not sensitive to the character of the external force since compressible motion creates shocks 
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and its divergent part decays faster than the solenoidal one [Ake Nordlund (2001), private 
communication], however, it may be sensitive to the presence of magnetic field that is known 
to help generate vorticity (Vazquez-Semadeni, Passot, & Pouquet 1996). 

These remarkable numerical observations lead us to a conjecture that the ratio 7 can 
be treated as a small parameter in the theory of 3D compressible turbulence. We assume 
that in the inertial region such turbulence is divergence-free, with the Kolmogorov time of 
velocity fluctuation, ti ~ l 3 ^ 2 , where / is the size of the fluctuation. Close to the dissipative 
range, shock structures start to play important role in energy transfer and dissipation. The 
turbulence in this region thus inherits certain properties of Burgers turbulence. The theory 
allowing to link the most singular, dissipative structures of turbulence with its velocity 
spectrum was suggested by She & Leveque (1994). This theory represents a turbulent 
cascade as an infinitely divisible log-Poisson process that has three input parameters. Two 
of these parameters are naive scaling exponents, 6 and A, of the velocity field and of the 
"eddy-turnover time", correspondingly: u\ ~ l @ , t\ ~ l A . The other parameter is the co- 
dimension, C, of the most singular dissipative structure (co-dimension is defined as dimension 
of space minus the dimension of the structure, C = d — D). The objective of the theory is 
to predict the so-called structure functions of the velocity field, defined as 

S p (l) = ([u(x + l)-u(x)] p )~lM, (2) 

where a is a component of the velocity field parallel or transverse to I. [According to 
the chosen component the structure functions are called either longitudinal or transversal. 
It is believed that both scale in the same way, so we do not specify what component is 
assumed in (2).] The velocity spectrum is a Fourier transform of the second-order structure 
function and is given by ~ k^ 1 ^^ . If the turbulent cascade depended only on local 
eddy interactions, then the naive Kolmogorov scaling of structure functions would hold, 
C(p) — p/3) an d we would recover the energy distribution E k ~ k~ 5 ^ 3 . Real turbulence is 
however intermittent, which means that its spectrum is not determined by the naive scaling. 
The She-Leveque theory predicts the scaling function as 

C(p) = e(i-A) P + c(i-s^ 3 ), (3) 

where £ = 1 — A/C. For the original derivation we refer the reader to the papers by She 
& Leveque (1994); She & Waymire (1995); Dubrulle (1994); more practical discussion can 
be found in (Grauer, Krug, & Marliani 1994; Politano & Pouquet 1995; Miiller & Biskamp 
2000). For 3D incompressible turbulence, the naive scaling exponents take the well-known 
Kolmogorov values = 1/3, and A = 2/3, and the dissipative structures are known to 
be filaments, so their co-dimension is C = 2. With these input parameters, formula (3) 
reproduces experimental results for incompressible Navier-Stokes turbulence with accuracy 
of several percent up to p = 10. 
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In our model of Kolmogorov-Burgers turbulence, the inertial range naive scaling ex- 
ponents are Kolmogorov ones, while the dissipative structures are quasi- ID shocks, which 
gives C — 1 and £ = 1/3. Formula (3) now reads 



This gives for the second-order structure function (uf) ~ /°' 74 , which reproduces the Larson 
law (Larson 1979, 1981), and the velocity power spectrum is given by E k ~ A; -1 - 74 , in 
a good agreement with numerical results by PN (2000). The intermittency correction to 
the Kolmogorov scaling is even larger for the first- order structure function, {\u{\) ~ l 0A2 , 
which can be checked observationally or numerically in an easier way. Our analysis here is 
analogous to the analysis of incompressible MHD turbulence by Grauer, Krug, & Marliani 
(1994); Politano & Pouquet (1995), and also by Miiller & Biskamp (2000) who noted that 
the most singular structures in such turbulence are micro-current sheets. The sheet-like 
dissipative structures together with the assumption that the energy cascade is given by the 
Kolmogorov rather than the Iroshnikov-Kraichnan mechanism, led the latter authors to the 
same prediction for the structure function scaling as our formula (4), which turned out to 
be in good agreement with numerical results. This indicates that both systems, though 
completely different, belong to the same class of universality, in agreement with the ideas 
put forward by Dubrulle (1994) and She & Waymire (1995). 



Our analysis relied considerably on sheet-like shock structures. Analogous considerations 
for the filament and core singularities would give E k ~ k~ imi for filaments (C = 2), and E k ~ 
^,-1.685 f or cores ((j — 3)_ w\ these spectra are steeper than the Kolmogorov one. Although 
the shock-like dissipative structures are clearly seen in simulations, the Mach number and 
the resolution are not large enough to make precise comparison with molecular clouds. As 
an important generalization of the theory, one can imagine that the dissipative structures 
are rather complicated on large scales and have dimensionality greater than two. This is 
consistent with observed fractal structures of the density distribution, whose dimensionality is 
close to D = 2.3 (Larson 1992; Elmegreen & Elmegreen 2001; Chappell & Scalo 2001). Since 
the substantial part of dissipation occurs in shocks, see, e.g., (Ostriker, Stone, & Gammie 
2000), the dimension of the most singular dissipative structures may be close to D = 2.3 as 
well. [It would be interesting to check numerically the degree of correlation of the density 
field and the dissipation field.] Substitution of C — 0.7 in our formula (3) leads to the 
first order structure function (\ui\) ~ /°' 55 , and to the energy spectrum E k ~ /c~ L83 . The 
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3. 



Generalizations 



- 6- 



results of recent observations, see, e.g. (Brunt & Heyer 2001), seem to agree well with these 
predictions. However, more precise measurement of the structure functions scaling would be 
required to indicate what structures are most important. The intermittency correction to 
a scaling exponent of the first-order structure function is large enough to be detectable in 
numerical experiments, but higher-order structure functions are more difficult to measure. 
An attempt to infer such structure functions from observations was made by Miesch, Scalo, 
& Bally (1999); Ossenkopf & Mac Low (2000), but the scaling was not established due to 
limited inertial ranges. 

Another important question is the relation of the obtained spectrum to the initial mass 
distribution function. The consideration of PN (2000) was based on an implicit assump- 
tion of the mean-field approximation, while our explanation of the observed steeper-than- 
Kolmogorov spectrum is essentially based on intermittency effects. In the presence of strong 
fluctuations, this relation may be modified, also acquiring intermittency corrections. 

I am very grateful to Richard Larson, Ake Nordlund, Annick Pouquet, Dmitri Uzdensky, 
and Enrique Vazquez- Semadeni for many valuable discussions, suggestions, and comments 
on both the physics and the style of the paper, and to Lars Bildsten, Bruce Elmegreen, 
Mordecai-Mark Mac Low, and John Scalo for useful remarks. 
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